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O ' We discuss the effects of the off-shell nucleon propagation for the particle produc- 

r'~'~) ' tion. Results from nonequilibrium Kadanoff-Baym evolution modeling a heavy ion 

collision are presented. The production rates are compared to equilibrium emission 
[^. t rates obtained in the same approximation for the self-energy. The comparison of 

semiclassical and quantum in-medium production rates is performed within the T- 
matrix approximation. Self-consistent T-matrix resummation allows us to discuss 
^. ■ also the effect of the off-shell propagation on the effective in-medium cross section. 
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^^ ■ 1 Nonequilibrium evolution 

o ; 

In medium modifications of nucleon properties are expected to be important 
at normal nuclear densities and higher. Most of these effects can be taken 
into account by a modification of the effective nucleon mass, cross sections 
and mean-field. In a strongly interacting medium one expects that off-shell 
propagation of nucleons is crucial. On the other hand, transport equations 

ri ' using on-shell quasi-particles proved to be very successful in the description of 

^ , intermediate energy heavy ion collisions. 

Daniclcwicz in the first work comparing an actual solution of the nonequi- 
librium Kadanoff-Baym equations with its quasi-particle limit found a slower 
equilibration in the quantum evolutional. Thus the relaxation rate obtained 
from the Boltzmann collision term is larger than obtained from the Kadanoff- 
Baym collision term, including off-shell propagation and memory effects. The 
modification of such global transport coefficients due to off-shell propaga- 
tion can be taken into account by an effective renormalizatian of the scat- 
tering amplitudes between quasi-particles. In a Fermi liquid □ this amounts 
to a multiplication of the cross section by the the factor z(pp ) 4 (z(p) = 
(1 — dlZeT,(p, ui)/duj\ LU=L j p )^ 1 ). For cold nuclear matter this leads to a re- 
duction of the cross section by a factor 2 or more. 

Production of very energetic mesons and photons has been observed in 
heavy ion reactions. This particle production is subthreshold, i.e. the energy 
available in single nucleon-nucleon collision is insufficient to produce such an 
energetic particle. Fermi motion of nucleon inside the colliding nuclei shifts 
the effective threshold to higher energies but cannot explain the high momen- 
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Figure 1: Production rate of particles of mass 140MeV. The solid and dashed-dotted lines 
represent the results with and without imaginary part of the contour in the meson self- 
energy calculation. The Fermion evolution is in both cases done with the imaginary part of 
the contour, i.e with the initial correlations. 



turn part of the spectra. The quasi-particle transport models invoke multiple 
collisions, or pion-nucleon collisions (with pions produced in earlier nucleon- 
nucleon collisions) to explain the production of energetic mesons or photons. 

Subthreshold particle production could be a testing ground for the off- 
shell effects in nuclear propagation. Off-shell spectral function of nucleons 
in hot nuclear matter allows for meson production in the simplest one-loop 
approximation for the meson self-energy a. Off-shell propagation of nucleons 
requires the solution of evolution equations for the two-time Green's func- 
tion, the Kadanoff-Baym equations. 16 years after the first calculations a it 
is still limited to spatially homogeneous systems. A numerical solution of a 
homogeneous transport equation could be used in the estimate of the particle 
production in heavy ion collisions because of the approximate scaling of the 
number of produced particles with the number of nucleons participating in the 
collision, even for the subthreshold production. The production rate in real 
collision could be estimated by scaling the production rate in a homogeneous 
system by a geometric. factor. 

In previous work 13 we studied the production in the nonequilibrium dy- 
namics of nucleons. The evolution of the system follows the Kadanoff-Baym 
equations, in the direct second Born approximation 

(i-^--u> p ^G<>\p,ti,t 2 ) = J dt'V(p,ti,t)G(p,t,ta) . (1) 



The integration contour C involves an imaginary part giving, within the second 
Born approximation, the ground state of nuclear matter with two boosted 
Fermi spheres. The real part of the contour gives the collision integral with 
memory. In order to get the spectral function far from the quasi-particle pole, 
small steps in times ti, £2 have to be usedB. 
n The meson production rate is obtained from the one-loop meson self-energy 

y 

where 



n<( P ,t 1 ,t 2 ) = -zA 2 y^G < (p-q,ii, 



h)G > ( P ,t 2 ,t 1 ) (3) 



and Dq is the vacuum meson propagator (we neglect the influence of the meson- 
fcrmion interaction on the nucleon evolution and we calculate only the produc- 
tion rate not the full meson transport equation). The results in Fig.l show 
that the production rate is oscillating with large amplitude at initial times. 
It also changes sign. It is a manifestation of the oscillating time exponent, 
which gives the energy conservation in the limit to ~^ —00. In an interacting 
system with off-shell nucleons the positive oscillations dominate giving a finite 
production rate from one-loop diagram. It is however not possible to define 
a particle production rate at initial times. The initial correlations take into 
account only the fermion-fermion interaction via a two-body potential. There 
is no initial meson content in the collision, i.e. the initial state is not correlated 
with respect to the meson-fermion coupling. This and the previously discussed 
meson formation time (oscillations of the rate) dominate the meson production 
at initial times. The meson production rate can be defined only at large times 
t ~ 20fm/c where the nucleon sector is basically equilibrated. 

2 Equilibrium rates 

Production rates in equilibrium finite temperature nuclear matter can be ob- 
tained with much less numerical effort and higher accuracy a. Equilibrium 
spectral functions for nucleons, at finite temperature and density, are obtained 
by iterative solution of the set of equations 

±»G<>](p,w) = A(p, «)/(«) [(1 - /(c))] , 

-2Tm£(p,w) 

A(p, w) = 



(u — p 2 /2m — 1ZeY>(p, lo)) 2 + ImS(p, u>) 2 
S<W(p,c)=^(G<,G>), 



7leT,(p, oj) = — 

J 7T 

d?pdio 



dw ImTi{p 1 uj 



(27T) 



G<( P ,L0)=p 



The functional T in the equation for S < [ > 1 corresponds to the direct second 
Born approximation as in the nonequilibrium evolution. In Fig. 2 we present 
the spectral function obtained at T = lOMeV. The spectral function is broad 
far from the Fermi momentum, but gets narrower for momenta close to it. As 
expected for a Fermi liquid the spectral function at the Fermi momentum gets 
very peaked when decreasing the temperature. 

The quantum production rate of mesons with momentum p and energy 
Clp = \J p 2 + M 2 in the one-loop approximation for the meson self-energy is 

S = 4A2 / Sf^ w)A{q ~ p,uj ~ np){i m]f{uj ~ Qp) - (4) 

The semiclassical production rate is 
dN(p) f d 3 px rf 3 p 4 



d?pdt 



/(pi) ... (1 - f( Pi )) , (5) 

with the matrix element for particle emission given by 

\M\ 2 = A\ 2 V{p 2 -p A f ( \G+(cj p1 - n„ Pl - q)\ 2 



+ \G + {oj ps +n q ,p 3 + q)\ 2 j 



(6) 



This matrix element is obtained by including self-energy corrections in the one- 
loop quasi-particle meson self-energy The sclf-cncrgy in the retarded propa- 
gators in \M\ 2 as well as the single particle energies u>i are taken from the 
quantum calculation of finite temperature nuclear matter. 

The production rate for particles of mass 140MeV is ~ 2 times larger in 
the quantum calculation. The self-consistent spectral function involves pro- 
cesses with multiple sequential nucleon-nucleon collisions which populate far 
off-shell regions of the spectral function. Also the one-loop self-energy with 
off-shell fcrmion propagators includes contribution from 3-nucleon collisions 
corresponding to self-energy insertions on both propagators in the loop. These 
contribution are not taken into account by the 2-body matrix element (Eq. 
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Figure 2: (left panel) Nucleoli spectral function at normal nuclear density and T = lOMeV 
for several momenta, (right panel) Production rate of mesons of mass 140 MeV at equilib- 
rium. Quantum and semiclassical production rates are represented by solid and dashed lines 
respectively. 



0). The difference between the semiclassical and quantum production rates 
increases with the energy of the produced particle. On the other hand, the 
emission of soft particles is stronger in the semiclassical calculation. All the 
emission rates increase strongly with temperature. Quantum production rates 
are similar as obtained from the noncquilibrium evolution U (Sec. [j]) at large 
times. The semiclassical rates published in the work El are different and obvi- 
ously wrong. The overall collision rate is larger in the semiclassical collision 
integral, similarly as in the noncquilibrium calculation. 



3 T-matrix approximation 

In this section we discuss in-medium modification of the cross-sections. This 
requires the use of the T-matrix approximation instead of the Born self-energy. 
T-matrixJadder resummation with off-shell propagators has only recently been 
achieved B. The equation for the T-matrix is 



<p|t ± (p,^)| p > = y(p, P ) 



d 3 k 



d 3 q 



V(p,k) 



(2tt) 3 7 (2tt) 3 

<k|g ± (P,^|q><q|r ± (P,c)|p > 

where the disconnected two-particle propagator is : 

< p|£±(P, fl)|p' >= (2vr) 3 ( 5 3 (p - p') / ^ / t(G<(P/2 +p,uj- u) 



(7) 
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Figure 3: (left panel) Single particle width obtained from Eq. (0) in the quasi-particle 
T-matrix approximation (solid line) and from the Boltzmann collision integral (dashed line), 
(right panel) Production rate of mesons as function of their mass in the quasi-particle T- 
matrix approximation. Quantum and semiclassical production rates are represented by solid 
and dashed lines respectively. 



G<(P/2-p,w')-G>(P/2 + p, 



LO — CO 



)G>(P/2-p,w))/(n-w±ie) (8) 



The T-matrix equation has to be solved iteratively together with Eqs. (0). 
The functional defining the self-energy is now 



< (p - k)/2|ImT+(p + k, n + w)|(p - k)/2 > A 



(9) 



instead of the Born approximation used in the previous section. The quasi- 
particle approximation in the T-matrix ladder consists in replacing the two 
fermion propagator (m by quasiparticle propagatorsB. First applications of the 
self-consistent T-matrix scheme to finite temperature nuclear matterUl show im- 
portant quantitative differences between the quasiparticle and self-consistent 
T-matrix results. Both at high temperatures where scattering is important 
and at lower temperatures close to the pairing transition. 

The quasi-particle T-matrix approximation gives in-medium cross sections, 
but all the nucleon propagators are on-shell in the ladder resummation. The 
single particle width and spectral function are the same as obtained using 
the Boltzmann collision integral with the cross section obtained from the in- 
medium T-matrix. Indeed the left panel in Fig. 3 shows the equivalence of 
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Figure 4: (left panel) Single particle width obtained from Eq. (H) in the self-consistent 
T-matrix approximation (solid line) and from the Boltzmann collision integral (dashed line), 
(right panel) Production rate of mesons as function of their mass in the self-consistent T- 
matrix approximation. Quantum and semiclassical production rates are represented by solid 
and dashed lines respectively. 



the two calculations (the small differences are due to numerical inaccuracy and 
approximate treatment, angular averaging, of the energy conserving S function 
in the collision term). The semiclassical production rate is defined as in the 
Born approximation except that the two-particle potential is replaced by the in 
medium T-matrix. The quantum production rate differs from the semiclassical 
one only by the inclusion of 3-nucleon processes corresponding to having both 
propagators in the meson self-energy loop off-shell. This contribution is most 
important at large meson energy and increases with temperature. 

The self-consistent T-matrix, using off-shell nucleon propagatorSjJn the 
ladder, leads to different spectral function and single particle energiesQ. In a 
self-consistent calculation there is no simple relation between the single-particle 
width obtained from the Boltzmann equation and from the T-matrix scheme. 
As in the Born approximation we find that on-shell propagation leads to larger 
relaxation rates. Part of this effect can be understood as due to the different 
densities of states for the two-particle propagators in the quasi-particle and self- 
consistent T-matrix calculations 3. The optical theorem, if using the density 
of states corresponding to on-shell nucleons, leads to an overestimation of the 
cross-section. 

The production rate is larger in the quantum calculation, especially for 
energetic particles. The difference comes from the 3-body contribution dis- 
cussed above and from the far off-shell regions of the spectral function. As in 
the Born approximation, the spectral function obtained in the self-consistent 
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iteration is larger (wider) far from the shell. On the other hand, close to the 
quasi-particle pole the spectral function calculated from the Boltzmann col- 
lision term is wider, reflecting the systematically larger single particle width 
obtained with on-shell nucleons. 

4 Summary 

We have calculated particle production in nonequilibrium quantum transport 
framework. The production rate could not be defined in the initial stage of the 
evolution. It is possible only at large times, when the system is already equili- 
brated. The corresponding rates can be calculated using equilibrium spectral 
functions, obtained in an iterative solution of coupled equations (|J). The 
quantum production rate is larger than the semiclassical one for energetic par- 
ticles. This can be understood as due to 3 and sequential n-nucleon collisions. 
In-medium T-matrix calculations allow to include in-medium modifications of 
the cross section in the picture. The cross section for Boltzmann type collision 
integrals and for the semiclassical production, rates has to be modified by a 
factor describing a different density of states 0. After such a renormalization 
of the cross sections the relaxation rates would be similar in the quantum and 
semiclassical transport equations, but the difference between the quantum and 
semiclassical production rates would increase. 
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